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Abstract. In this article we present some recent results on identifying correctly the relativistic 
multipole moments of numerically constructed spacetimes, and the consequences that this 
, correction has on searching for appropriate analytic spacetimes that can approximate well 

the previously mentioned numerical spacetimes. We also present expressions that give the 
quadrupole and the spin octupole as functions of the spin parameter of a neutron star for 
various equations of state and in a range of masses for every equation of state used. These 
results are relevant for describing the exterior spacetime of rotating neutron stars that are made 
up of matter obeying realistic equations of state. 



1. Introduction 

The multipolar expansion of the gravitational field in Newtonian gravity is a straightforward 
exercise. In spherical coordinates one could express the field as an expansion in powers of the 
radial coordinate and the appropriate angular eigenfunctions of the Laplace operator. One could 
be tempted to straightforwardly apply the same procedure to General Relativity (GR), but in 
this case things are more complicated due to the non-linear nature of the theory (the principle 
of superposition does not apply in GR). The task of defining the multipole moments in GR 
was undertaken in the beginning of the 70s by Geroch and Hansen who defined the multipole 
moments of an asymptotically flat spacetime in the static and stationary case as tensors at 
infinity [HE]- Alternative definitions of the relativistic multipole moments were also given in 
the early 80s by Simon [3] and by Thorne [3] which were closer to the spirit of the Newtonian 
asymptotic expansion, where Simon's moments end up to be the same moments as Geroch's and 
Hansen's, while Thome's moments are coordinate dependent and an appropriate choice of the 
coordinates should be made in order to read the usual Geroch-Hansen moments (for a review 
see [S]). Finally in 1989, Fodor et al. [BJ developed an algorithm for computing the multipole 
moments of a spacetime that is additionally axially symmetric, taking advantage of the insightful 
Ernst-potential formalism, expressing thus the scalar (in this case) multipole moments in terms 
of the Ernst potential. 

In the case of stationary and axially symmetric spacetimes, the relativistic multipole moments 
can fully characterize the spacetime as one can see in [BJ. Thus the multipole moments are of 
interest as quantities that capture all the properties of a given spacetime. On the other hand, 
they are also of interest on experimental grounds, since they can be measured in gravitational 



wave signals from gravitational inspirals. In particular, Ryan (7| related various observed 
properties of these inspirals to the multipole moments of the background spacetime. Specifically 
he gave relations that connected the multipole moments of a spacetime to precession frequencies 
of orbits on that spacetime that deviate slightly from being circular equatorial ones as well as 
to the number of cycles of the gravitational wave signal that is emitted by a low mass object 
inspiraling adiabatically into such a spacetime. These same expressions can be also used in 
the context of accretion discs and the quasi-periodic oscillations (QPOs) observed in the X- 
ray spectrum of low mass X-ray binaries (LMXRBs) [8], where one, assuming a relativistic 
precession model for the QPOs [S], can relate the QPO frequencies to the multipole moments 
of the central object [10J. The connection of the multipole moments to the structure and the 
physical properties of a compact object, like a neutron star, was first attempted by Laarakkers 
and Poisson, who related the multipole moments of a neutron star to the equation of state (EOS) 
for the matter in its interior 

Furthermore, the multipole moments can be used to construct analytic spacetimes that will 
have prescribed properties. Thus, in the case that we would like to have an analytic description 
of the spacetime around a neutron star with specific physical properties, we could use the 
appropriate multipole moments to construct such a spacetime following the procedure presented 
in [12\ [T3l [T^l [T5l [T6] (the multipole moments can also be important in different approaches to 
the problem of constructing geometries for neutron stars, such as [TTJ [18]). 

Working on the previously mentioned applications of multipole moments to analytic and 
numerical spacetimes, we came to realize that there is a problem with the way the multipole 
moments of numerically constructed spacetimes are identified from their asymptotic behavior. 
These issues were addressed and clarified in [19] where the correct expressions for the quadrupole 
and the spin octupole were derived using Ryan's [7] algorithm for relating the gravitational-wave 
spectrum AE (the energy emitted per unit logarithmic frequency interval) of a test particle 
that is orbiting on a circular equatorial orbit in an asymptotically flat, stationary and axially 
symmetric spacetime, to the multipole moments of that spacetime. In [19] various consequences 
of correcting the calculation of the numerical multipole moments were also discussed, such as 
the effect this correction has on attempting to approximate the numerical spacetime exterior to 
neutron stars. In particular it was shown that when one uses an analytic spacetime that depends 
on a number of parameters that can be connected to the multipole moments of a numerical 
spacetime, such as the analytic spacetime of Manko et al. [20^ 121]. in order to approximate that 
numerical spacetime, the use of the corrected values for the multipole moments improved the 
performance of the analytic spacetime and resulted to a better fit to the numerical spacetime. 

It was also shown, extending a previous result by Laarakkers and Poisson [11] , that the 
reduced quadrupole, q = M2/M 3 , and the reduced spin octupole, S3 = S3/M 4 , follow 
respectively quadratic and cubic dependance on the spin parameter, j = J/M 2 , just as in 
the case of the Kerr spacetime where these moments are — j 2 and — j 3 respectively, with the 
difference being that the proportionality constant is larger than 1 for neutron stars. 

In this work we briefly present these results, i.e., we give a short derivation of the correct 
expressions for the multipole moments and then present the resulting improvement in the fit of 
the numerical metric by the 3-parameter analytic metric of Manko et al. In addition we present 
some results from approximating the numerical metric with the 4-parameter two-soliton metric, 
that is extensively discussed in [16], that further strengthen our arguments in favor of using 
the corrected multipole moments as criteria for matching an analytic spacetime to a numerical 
one. Finally we give the proportionality constants between the quadrupole and the square of the 
spin parameter, j 2 , and the spin octupole and the cube of the spin parameter, j' 3 , for different 
equations of state and for different masses. These coefficients are directly comparable to the 
coefficients of Laarakkers and Poisson for the quadrupole. All the numerical spacetimes are 
calculated using the RNS numerical code [22J. All the expressions are in units of G = c = 1. 



2. Asymptotic expansion of a metric and the multipole moments 

The line element used for the spacetime for the interior and the exterior of a neutron star, 
which has symmetry with respect to time translations and rotations, i.e., it is stationary and 
axisymmetric, is usually written in quasi-isotropic coordinates in the form, 

ds 2 = _ e 2» dt 2 + r 2 gin 2 6B 2 e -2u {d(j) _ udt) 2 + ^-u) ^2 + (1) 

where v, B, u, and C, are the four metric functions, all functions of the quasi-isotropic coordinates 
(r,9), as it was introduced by Butterworth and Ipser [23]. From Einstein's field equations, one 
can show that the asymptotic expansion of the metric functions takes the following form, 
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In the formulae above Pi are the Legendre polynomials expressed as functions of /j = cos 9, T t 
are the so called Gegenbauer polynomials (similar to the Legendre polynomials, also functions 
of fi), and M, J are the first two multipole moments (the mass and the spin) of the spacetime. 
The rest of the coefficients are related to the higher multipole moments. 
These metric functions can be rewritten using the following redefinitions, 
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Then the combinations of v and /3, 
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along with uj could again be expressed as power series in 1/r (this is done in [241 1251 [22]) i n t ne 
same manner as in Eqs. (|2l3l4p : 
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By a simple comparison between the above expansion and the corresponding ones in Eqs. (|2l3l4p . 
one can see for example that M| = — z>2 and S3 = 1^2. The coefficients M 2n and S 2n -i were 
mistakenly identified as the mass and current-mass moments, respectively, of the corresponding 
spacetime and it is exactly these quantities that the numerical code of Stergioulas and Friedman 
[22], RNS, provides. 

In order to identify the correct multipole moments, i.e., the Geroch-Hansen multipole 
moments, of a spacetime given by the line element (pQ), we have used Ryan's expression [7] that 
relates the multipole moments with the energy change per logarithmic interval of the rotational 



frequency AE = —dE/dlogQ for circular equatorial orbits in a stationary and axisymmetric 
spacetime. Thus, for the line element of Eq. ([1]), we first expressed the orbital frequency of 
circular equatorial orbits, 
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as a power series in x = (M/r) 1 / 2 and then inverted it to obtain, 

v 3 jv A 1 5 jv 6 (97 + 286 + 56 j 2 - lUq)v 7 



v + — + J — + — (13 + 46 - 6q)v b + — + 



2 3 24 v 2 144 

(373j + 2926j - 330jg - 270w 2 )v 8 9 

360 [V h 1 ' 

where v = (Mil) 1 / 3 , j = q = jjfc, w 2 = jfc, b = Jj£. Then we expanded 
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which is the energy per unit mass for a specific circular orbit, as a power series in x = (M/r) 1 ^ 2 
and substituted the previous expression for x. From the resulted expression we calculated AE, 
from the equation, 

A p_ d E _ vdE h a\ 

and arrived to the expression, 
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If one compares the previous expression to the one produced by Ryan, 



A ~ 1 2 1 4 20 Si s / 27 M 2 \ 6 28 Si 7 

AE = -v 2 v A H ^y 5 + 1 % )v b H W 

3 2 9 M 2 V8 M 3 y 3 M 2 

/ 225 80 S 2 70 M 2 \ 8 / 81 Si SxM 2 S 3 \ „ 
+ ("16" + 27 M^ + ¥m^T + Ui^ + 6 l^- 6 iiFj t; + - (15) 

where Si = J, one can see that from the coefficients of v e and v 9 terms of the two series, the 
following values for the quadrupole and the spin octupole can be obtained: 

M? H = -z> 2 -^Q + 6)m 3 = M 2 *-^Q + 6)m 3 , (16) 

S$ H = fa* Q + ^-M 4 = S 3 *-f (\ + b)jM\ (17) 



Henceforth we will omit the superscript , that indicates the Geroch-Hansen moments, in 
M 2 , S3 when we refer to the correct multipole moments. 



Laarakkers and Poisson performed in [TT] a similar calculation in order to identify the 
quadrupole of the metric ([1]), but their result was missing the last term in Eq. (|16p . i.e., the 
term with (1/4 + b). That was because they erroneously assumed that the correct asymptotic 
behavior for the metric should be that of Schwarzschild, which has B = 1 — M 2 /4r 2 , which 
corresponds to b = —1/4. Thus they had a priori fixed the correcting last term to be zero. 
Under a more careful consideration though, there is no reason to fix the asymptotic behavior of 
a stationary and axially symmetric spacetime in this way. A counter example for this is the case 
of the Kerr spacetime, where one can see that in quasi-isotropic coordinates the metric function 
B is given by, 

B KclT - 1 -p , (18) 



which corresponds to having 6 Kcrr = — (1/4)(1 — j 
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3. Improvement in approximating a numerical spacetime with analytic spacetimes 
due to the correction in the moments 

We will now present the effect that the correction in the moments has when one attempts to 
use an analytic spacetime to approximate the exterior space time of a neutron star, that is 
calculated numerically, by matching the multipole moments of the analytic spacetime to the 
multipole moments of the numerical spacetime. Berti and Stergioulas [26J tried to match a 
three-parameter analytic solution, the solution of Manko et al. [20], to a wide diversity of 
uniformly rotating neutron-star models. Each analytic solution was constructed so that its first 
three multipole moments were equal to the corresponding moments (mass, spin and quadrupole) 
of the particular numerical neutron star, where these moments were read from the corresponding 
numerical metric. For their calculations, they used as moments the quantities M% n and S^n-i- 
Their conclusion was that this type of analytic solution was quite good to describe the external 
metric of all kinds of fast rotating neutron stars. Since the specific metric cannot assume low 
values of quadrupole moment, the metric is not adequate to describe rotating neutron stars with 
rotation lower than some value. We have used the Manko et al. solution to demonstrate the 
effect of correcting the multipole moments, compared to the results by Berti and Stergioulas. 

Since the Manko et al. solution is in the previously mentioned sense handicapped, and since 
it has only three free parameters which provide the freedom to explore the behavior of the 
metric by changing only the quadrupole, we have also used another analytic solution, the so 
called two-soliton solution [27], that can also be used to approximate the exterior spacetime of 
neutron stars [16J and has four parameters, which allow into play the spin octupole as well. Thus, 
with the freedom now to vary both the quadrupole and the spin octupole, we have performed 
comparisons between the two-soliton and specific numerical metrics in order to see which values 
of multipole moments give the best fit. 

In both cases, for the Manko et al. and the two-soliton, we have found that when the 
parameters of the analytic metric are calculated using the correct values for the multipole 
moments, the fitting of the analytic metric to the corresponding numerical one is better. This 
result also supports our claim that the way that an analytic metric should be matched to the 
numerical metric is by identifying the corresponding multipole moments of the two spacetimes. 

In Figure [T] we present a typical comparison between the analytic metrics that are calculated 
using the correct moments, i.e. M2, S3, and the analytic metrics calculated using the previously 
assumed moments. Specifically the figures show the relative difference between the numerical 
and the analytic metric ((gfj — gf^/gfj) as a function of coordinate radius for the Manko et al. 
analytic metric and the two-soliton analytic metric, calculated for the values M|, S3 and for the 
values M2 , 53 (the value of the spin octupole is relevant only for the two-soliton spacetime which 
has 4 parameters). The red dashed curves correspond to the Manko et al. metric produced from 
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Figure 1. A typical log- log plot of the relative difference between the numerical and the analytic 
metric ((g™- — gfj)/gij) for a specific numerical model (model #19 of EOS FPS of [E]). The left 
plot is for g tt and the right one for g t §. The red dashed curves correspond to the Manko et al. 
metric produced from M| and the blue dashed curves correspond to the Manko et al. produced 
from Mi. The orange solid curves correspond to the two-soliton produced using the moments 
M|, S3 and the green solid curves correspond to the two-soliton produced using the corrected 
moments M2 , -S3 . 



M| and the blue dashed curves correspond to the Manko et al. produced from M%. The orange 
solid curves correspond to the two-soliton produced using the moments M^S^ and the green 
solid curves correspond to the two-soliton produced using the corrected moments M2,Ss- The 
left plot in Figure Q] shows the relative difference for the gu metric component while the right 
plot shows the relative difference for the gts metric components. 

In order to further test the effect that varying the moments M2 and S3 has on the performance 
of the analytic metric to approximate the numerical metric, we have used the two-soliton metric 
to investigate the parameter space of the variation of the moments. As a measure of the ability 
of the analytic metric to approximate the numerical metric we have used the quantity, that we 
call, "overall mismatch" between the analytic and the numerical metric functions, that is defined 
(see [19]) as 




where Rs is the radius r at the surface of the star. For each numerical model of the various 
models of uniformly rotating neutron stars that we have constructed, we formed a set of two- 
soliton spacetimes that have the same mass M and angular momentum J with the numerical 
model, but the quadrupoles and the current octupoles take the values M2 = M2(l — (IM2) and 

S 3 = 53(1 — dSs) respectively, where the (IM2 and dS% take various values. The quantities dM2 
and dSz denote the fractional differences of the corresponding moments used in the calculation 
of each two-soliton spacetime from the moments of the numerical spacetime. For each one of 
the two-soliton spacetimes belonging to the previously mentioned set, we calculated the "overall 
mismatch" between the analytic and the numerical metric functions producing contour plots of 
the "overall mismatch" on the space of dM2 and dSa for the different numerical models (one can 
find further discussion on this in [16J). Examples of these contour plots are shown in Figure [2j 
These contour plots show that the choice of multipole moments for an analytic spacetime 
that give the best approximation to the numerical spacetime, are those that correspond to the 
moments given by the equations (|16|17p . i.e., those with (dM2,dSs) ~ (0,0). 
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Figure 2. Contour plots that demonstrate that the appropriate moments to be used for 
matching an analytic spacetime to a numerical spacetime are the corrected moments il^S^. 
The plots show the contour plots of the overall mismatch between the analytic and the 
numerical metric for the tt (red curves) and tcj> components (blue curves) as a function of the 
fractional deviation of the quadrupole, (IM2, and the current octupole, dSs, of the analytic metric 
from those calculated from the numerical metric as they are given by the corrected expressions 
(assuming the same mass and angular momentum). The value of the overall mismatch increases 
from the inner to the outer contours. Because the contours of the crtt (almost horizontal) are 
orthogonal to the ones of er^ (almost vertical), the combination indicates an optimum choice for 
the multipole moments of the analytic space-time. That choice is the moments that have zero 
deviation from the moments of the numerical space-time. These contour plots correspond, from 
left to right, to the models #22 and #28 for the AU EOS and #28 for the L EOS of the models 
presented in [19] . 

4. Relation between the multipole moments and the spin parameter j 

The connection of the higher moments of a compact object with its spin parameter j was 
first attempted by Laarakkers and Poisson [11]. Specifically what they showed was that the 
reduced quadrupole moment of neutron stars constructed using realistic equations of state was 
proportional to the square of the spin parameter and the proportionality constant depended on 
the mass of the neutron star and the equation of state, i.e., 

q ~ -a(M, EOS)j 2 . (20) 

In |19] it was shown that this relation is true, even after the correction of the quadrupole 
moment because the correcting factor (1/4 + b) is proportional to j 2 , which is strikingly similar 
to the way that this factor behaves in the case of the Kerr spacetime. A similar result was shown 
to apply for the spin octupole as well, i.e., the spin octupole has a dependence on the cube of 
the spin parameter of the form, 

s 3 ~-(3(M,EOS)j 3 . (21) 

The calculations performed in [19) . were performed using evolutionary sequences, so these 
results were not directly comparable to the results of Laarakkers and Poisson. Here we perform 
these fits for the reduced quadrupole and the reduced spin octupole as functions of j 2 and j 3 
respectively and calculate the coefficients a and (3 for neutron star models of constant mass and 
varying rotation. The results are presented in Table [TJ 

These results are comparable to those presented by Laarakkers and Poisson in Table VII of 
[TT) . So, we can compare the quadrupole coefficients for the equations of state FPS and L that 
are given both here and in [TT]. We should first note that there is a typo in Table VII that 



Table 1. The table gives the coefficients a, P of q = —aj 2 and S3 = — /3j 3 for various equations 
of state ranging from soft (on the left) to stiff (on the right) and for different masses in the range 
of 0.9 — 2.1M . More details on the equations of state used can be found in [251 ES] 
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shows the fit parameters. The value of a for the 1.4M models for EOS FPS should be 4.2. 
To compare the results in our Table [1] and the results in Table VII of [TTJ we will calculate the 
relative difference between the quadrupole moments one would calculate from our coefficients 
and from the coefficients in [TT]. The relative difference in the quadrupole moments would be 

AM2 = M 2 = (-aj*M3) = ~T~ ' (22) 

i.e., equal to the relative difference in the coefficients. The results are shown in Table where 
we see that there is an increase in the error of calculating the quadrupole as we go to models of 
higher mass, that for the FPS EOS can be as high as 11% for the models of 1.8Mq. 



Table 2. The table gives the relative differences in the coefficients a for the quadrupole between 
the results presented in Tableland in Table VII of [11] for the equations of state FPS and L. 
The relative differences are given as %. 
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5. Conclusions 

In this work we have briefly presented some recent results (see |19| ) regarding the correct 
identification of the multipole moments of a spacetime given in quasi-isotropic coordinates, 



with application to the evaluation of the multipole moments of numerical spacetimes. We 
have also briefly discussed the effect that this correction has on using the multipole moments 
as parameters for constructing parameterized analytic spacetimes that approximate numerical 
spacetimes around compact objects. Specifically we showed that when the correct moments 
are taken into account, an analytic spacetime gives a better approximation of the numerical 
spacetime and further more, for the case of the two-soliton spacetime, the best approximation is 
the one that implements the corrected quadrupole and spin octupole. Finally we have discussed 
the fact that the quadrupole and the spin octupole seem to be proportional to the square of the 
spin parameter, j 2 , and the cube of the spin parameter, j 3 , respectively, and we have presented 
the proportionality constants for various equations of state and various masses. 

We hope that the coefficients presented in Table [TJ will be useful to anyone that would like to 
construct parameterized analytic spacetimes for the exterior of neutron stars, that involve the 
first four non-zero multipole moments or relations between them as parameters. 
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